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Abstract. 

Structural and dynamical fingerprints of evolutionary optimization in biological 
networks are still unclear. We here analyze the dynamics of genetic regulatory networks 
responsible for the regulation of cell-cycle and cell differentiation in three organism or 
cell types each, and show that they follow a version of Hebb’s rule which we term as 
coherence. More precisely, we find that simultaneously expressed genes with a common 
target are less likely to conflict at the attractors of the regulatory dynamics. We then 
investigate the dependence of coherence on structural parameters, such as the mean 
number of inputs per node and the activatory/repressory interaction ratio, as well as on 
dynamically determined quantities, such as the basin size and the number of expressed 
genes. 
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List of abbreviations 

GRN Gene regulatory network 

Th T-helper lymphocyte 

1. Introduction 

Gene regulatory networks (GRNs) constitute the backbone of intracellular functional 
organization at the molecular scale. These interaction networks are key to understanding 
the clockwork operation of a cell’s life cycle [T], mechanisms of response to environmental 
changes [2], robustness against random fluctuations pm, the effects of mutations pm, 
embryonic development in higher organisms [2j, etc. Thanks to an enormous amount of 
data generated by recent experimental and computational efforts, we now have access 
to gene expression profiles in continuous time and can use them to deduce underlying 
regulatory interactions pm, even speculate on the evolution of such interactions in 
historical time scales mm- 

Inferring the global GRN of an organism from time-resolved gene expression 
data is an ongoing challange [12|, H8] , Therefore, past decade witnessed a growing 
interest in identifying key principles that govern the structural organization of the 
GRNs US m IS, 17-1- Viewing the regulatory network as a collection of functional 
subunits is a popular paradigm m m USE ED G2], supported by the observation 
that certain motifs are frequently encountered in the regulatory networks of many 
organisms [23]. The GRN structure is ultimately determined and constrained by the 
requirement that the regulatory dynamics delivers a timely production of necessary 
proteins. The fact that some of the frequently encountered motifs promote dynamical 
stability and robustness to minor failures is therefore not surprising m ■ Controllability 
recently emerged as another defining feature of these complex systems [2.5], stressing the 
requirement for a better understanding of the interplay between the network architecture 
and its dynamical behavior. Despite the success of such approaches, it has been 
pointed out that there is need to develop new methods taking different edge signs into 
account [26f. Present investigation of coherent regulation in biological networks is a 
progress in this direction. 

Coherent regulation: Protein production constitutes about one-half of raw material and 
energy consumption within a growing bacterial cell and one-third for a differentiating 
mammalian cell [2ZU2HU2S]. Therefore, it is plausible to ask whether the gene regulation 
hardware is wired in a way to achieve the desired functionality with minimal use of these 
resources. Considering the structure-dynamics relation from the perspective of energy 
efficiency, we propose and provide evidence that the GRN architecture has been partly 
shaped to promote unity of purpose among simultaneously expressed genes sharing a 
common regulatory target. We refer to such cooperative action of regulatory genes as 
“coherent regulation” [30J. 
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The idea that the evolutionary pressure for economy may have shaped regulatory 
interactions is not new; for example, it has been exploited earlier to identify the class 
of Boolean functions that better model regulatory dynamics ra, or to investigate the 
frequency of gene duplication in microbes [32]. We claim that, network structures with 
energy-optimal functionality should be wired to suppress the expression of “opposing 
minority” regulators. These are transcription factors which, if expressed, would oppose 
but not significantly alter the target gene’s fate due to outweighing regulatory pressure 
favoring the status quo. Networks where such minority influences are suppressed 
would display a disproportionate degree of consensus among simultaneously expressed 
regulatory elements acting on a common target, i.e., exhibit coherent regulation. 

Note that, the definition of coherent regulation here is different from that used in 
the context of robustness analysis [33]. Yet another use of similar terminology appears 
in the categorization of network motifs [34], where the coherence of a motif is determined 
according to the compatibility of alternative directed paths connecting two nodes. In 
contrast, the degree of coherence defined in the present work is not only a function of the 
interactions (edges) in the network, but also of the expression states of genes (nodes). 

Investigation of coherence on biological networks requires information about the 
regulatory machinery in the cell; in particular the architecture (say, in the form of a 
directed graph) and the character (activation/inhibition) of interactions, as well as a 
detailed knowledge of the regulatory dynamics. Dynamical aspects of genetic regulation 
have been investigated both on small motifs composed of a few genes [35], and on larger 
networks [36,03]. Depending on the desired resolution, Boolean models [38, 391 [401 41], 
Petri nets [i2i m m, and differential equation based continuum models na M\ are 
the typical approaches employed for this purpose. A continuum model is indispensable 
for a high (time-)resolution study of regulatory dynamics. Simpler Boolean models 
have also found a wide area of applicability, mostly in studies where a coarse 
characterization of the (quasi-)static stationary states is acceptable [47, [39]. These 
approaches have been successful in modelling the regulation of cell cycle [38, SB, 49], 
cell differentiation [[50, 51, 52], circadian clocks [53] 54], etc. We test our hypothesis on 
Boolean systems due to their simplicity and accessibility, although our approach can be 
generalized in a straightforward manner to continuum models. 

The organization of the paper is as follows: Section [2] is a formal introduction to 
the Boolean network dynamics and coherent regulation. In Section 01 we introduce 
six regulatory networks of different organisms or cell types, for which well-established 
Boolean models of regulation were adopted from the literature. Section 0] reports our 
results which suggest a bias towards high coherence in these systems, upon comparison 
with appropriately constructed random networks. Section 0] investigates structural and 
dynamical features associated with coherent regulation. Finally in Section El we discuss 
our findings and provide motivation for further investigation of coherence in complex 
networks. Overall, the current work extends our earlier observation on a single GRN Pi 
to multiple organisms or cell types and suggest that a bias towards high coherence may 
be a generic feature of gene regulation in biological systems. 
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2 . Computational framework 

2.1. Time evolution 

Following the standart notation, we describe the expression level of a gene by a binary 
variable with values 0 (silent) or 1 (expressed). Therefore, the state of a GRN composed 
of n nodes at a given time t is a binary vector cr(t) of length n, where ajt) is the 
state of the i th gene (1 < i < n) at time t. Time evolution in a deterministic setting 
(which applies to all the models considered here) is described by an evolution operator 
T:cr(t + 1) = T[<x(f)]. Given and initial condition cr(0), the ultimate fate (steady state) 
of the system is a cycle of length q, where a member state cr* of the cycle satisfies the 
condition cr* = T q [cr*]. A fixed point is a trivial cycle with q = 1. 

The rules of dynamics for the biological networks considered in this study were 
adopted from the respective references [ill ESI EMESES El!- Therefore, the evolution 
operator T is different for each network, except the cell-cycle models of the two variants 
of yeast for which the same majority rule was employed. We provide network-specific 
details of the GRN structure/dynamics in Section El 

2.2. Quantifying coherent regulation 

For convenience, let us describe the regulatory influence of a node i on node j by 
Cij G { — 1, 0,1}, where ±1 indicates positive/negative regulation and 0 indicates absence 
of either. In order to assign a coherence coefficient a c to a GRN, we first define a single 
node i of the network to be coherently regulated in state cr if and only if the inputs 
from its “on” neighbors j (oj = 1) are all activatory or all repressory (see Fig. Q]), i.e., if 



(1) 


j 


j 


The degree of coherence for a state cr can then be measured by the fraction of coherent 



Figure 1: An example network with coherently and incoherently regulated genes. Black 
and white circles indicate active and inactive nodes, respectively. Node 4 receives 
coherent input (only repressive) since node 1 is inactive. But node 1 itself is incoherently 
regulated since it receives conflicting inputs from nodes 2&3. 
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nodes in that state: 
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A node with no input (both the numerator and the denominator vanish above) is defined 
to be perfectly coherent. 

We quantify the coherence of the whole network through its steady states 
(attractors) corresponding to the fixed points or limit cycles of the dynamics, in which a 
cell spends most of its time. Labelling different steady states with s, the global coherence 
coefficient of the network is expressed as 


a. 


= E 


9s 

w s r (s)i 

— E“c |< t . 1 

HS i = 1 


( 3 ) 


where q s is the length of the limit cycle s and w s is a weight (such as the relative basin 
size, i.e., the fraction of initial states that end up in the given attractor) subject to the 
condition '^f /s w s = 1. 


Null-model ensembles: Bias towards coherent regulation in a GRN can be assessed by 
comparing a c for the given network with the distribution of the same quantity in a 
representative ensemble of similar networks. We construct such a reference ensemble 
separately for each regulatory network detailed in the next section. The ensemble 
networks were generated by shuffling the edges of the original GRN sufficiently many 
times such that, the resulting network 

• is a connected graph, 

• strictly conserves the self edges along with the number of incoming and outgoing 
activatory/inhibitory interactons separately for each node , 

• statistically has no correlation with the original network, except for the local 
similarities imposed by the two constraints above. 

By construction, each gene in these random ensembles is locally subject to the 
same number of repressory and activatory inputs as in the original network, albeit from 
possibly different regulatory partners. Under the same rules of dynamics the fixed points 
and the corresponding a c values are generally different for each random network, since 
they are determined by the new global structure. For each ensemble, we generated 
10 4 non-isomorphic networks and calculated their coherence coefficient both with and 
without a basin-size dependent weight (w s ) assigned to each dynamical attractor. 


3. Investigated genetic regulatory networks 

In this paper, we investigate the degree of coherent regulation in six GRNs associated 
with different organisms or cell types: cell-cycle networks in Saccharomyces cerevisiae 
(budding yeast) [41], Scizosaccharomyces pombe (fission yeast) [ 38 ] and mammals [ 32 ], 
cell-differentiation networks of Arabidopsis thaliana whorls [50], Th lymphocyte [58] and 
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(a) 



(d) 



Figure 2: Gene regulatory networks for (a) Fission yeast cell cycle [2Sj; (b) 
Mammalian cell cycle [HH]; (c) Budding yeast cell cycle m, (d) Arabidopsis 
thaliana whorl differentiation 133; (e) Myeloid differentiation [HE]; (f) Th-Lymphocyte 
differentiation [5BJ. Here green and red edges indicate activation and inhibition 
respectively which arise through physical or genetic interactions between the nodes 
representing genes, proteins, protein complexes or transcription factors. 


myeloid progenitors [56]. A graphical description of each GRN is given in Fig|2l These 
dynamical models were chosen from the literature, subject to the criterion that they re¬ 
produce the experimentally observed steady-state expression profile (s) after truncation 
to Boolean variables. Below, we give a brief description of each GRN. 

Scizosaccharomyces pombe (fission yeast) cell cycle: The fission yeast cell-cycle net¬ 
work was modeled by Davidich and Bornholdt [3H] as a network with 10 nodes (Fig. EK)- 
The dynamics is governed by threshold functions (see TableED which yield 12 fixed points 
and a fixed cycle. The fixed point with the largest basin matches the biological G i phase. 

Mammalian cell cycle: Faure et al. [29] analysed the regulation dynamics with syn¬ 
chronous, asynchronous, and mixed updating schemes in this GRN model composed of 
10 key regulatory elements. Regulation dynamics is given in terms of logical expressions 
(see Table Hj) which are determined according to available experimental evidence for each 
node. The resulting dynamical attractors are independent of the updating scheme and 
include a fixed point and a limit cycle, in agreement with the experimental expression 
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data. Note that, the visual depiction of the GRN given in Ref.[39] is inconsistent with 
the used logic update functions. We here remained faithful to the given logical expres¬ 
sions, after verifying that they reproduce the reported steady states. The structure of 
the GRN consistent with the interactions in Table 0 is given in Fig. 0). 

Saccharomyces cerevisiae (budding yeast) cell cycle: The model proposed by Li et 
al |H! is composed of 11 nodes (FigJ2b). This popular model reproduces the Gi phase 
as the dominant attractor of a simple Boolean dynamics, as well as the intermediate 
phases of cell division. Time evolution is governed by threshold functions given in Ta¬ 
ble El which yield 7 fixed point attractors. The attractors other than G\ have relatively 
small basins and, to our knowledge, no clear biological interpretation. 

Arabidopsis thaliana whorl differentiation: Mendoza et al. |55] use the network in 
Fig. 01 in order to model the dynamics of flower morphogenesis. The interactions in 
the 11-node model network are again inferred from experimental data. Different initial 
conditions evolved by the proposed rules of dynamics yield 6 point attractors, 4 of which 
have a clear biological interpretation. 

Myeloid differentiation: A Boolean model was set up by Krumsiek et al. [56] in order 
to understand the mechanisms underlying myeloid differentiation from common myeloid 
progenitors to megakaryocytes, erythrocytes, granulocytes and monocytes. The 11-node 
model network (Fig. [2b) is composed of relevant transcription factors which evolve (in 
time) under separate logical update functions, again inferred from the available experi¬ 
mental evidence. The dynamics gives rise to 5 point attractors, where 4 are in agreement 
with microarray expression profiles of the mature cell types. It is pointed out that the 
fifth attractor cannot be realized during physiological hematopoietic differentiation. 

Th-Lymphocyte differentiation: Remy et al. m proposed this regulatory network 
model for the differentiation of T-helper lymphocytes (ThO) cells into Thl and Th2 in 
the vertebrate immune system. Discrete-time evolution of each node in the model net¬ 
work is governed by node-specific logical rules given in Table 0 The network structure 
in Fig. 0 is dominated by activatory interactions. The dynamics settles into 3 steady 
states, in agreeement with the gene expression levels in the ThO, Thl and Th2 cells. 

It was shown earlier that structural parameters such as the mean number of 
incoming edges per node (k) and the fraction of up-regulating interactions (p) in the 
network are important determinants of global coherence, while the size of the network for 
fixed (k,p) only weakly influences a c |3()j . Observed values of (n, k,p ) for each network 
are listed in Table 0 For the GRNs we consider here, k varies within [1.67,3.1], while 
p changes in the interval [0.32,0.65] (self-edges are excluded). All models have roughly 
the same network size (~ 11 nodes). 
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Organism/Cell type 

Func. 

Network parameters 

# of attr. 



n 

k 

P 

all 

bio. 

S. pombe 

C.C. 

10 

2.3 

0.35 

13 

1 

Mammalian cell 

c.c. 

10 

3.1 

0.32 

2 

2 

S. cerevisiae 

c.c. 

11 

2.64 

0.51 

6 

1 

A. thaliana whorls 

diff. 

11 

2 

0.55 

6 

4 

Myeloid progenitor 

diff. 

11 

2.36 

0.42 

5 

4 

T h- ly mpho cy t e 

diff. 

12 

1.67 

0.65 

6 

2 


Table 1: Structural/dynamical features associated with the six regulatory networks 
adopted from the literature, “c.c” and “diff.” refer to regulation of cell cycle and 
differentiation respectively. Last two columns exclude the trivial point attractor cr = 0. 


4. Coherence in biological networks 

Table El and Fig. El summarize the outcome of our coherence analysis on the GRNs 
depicted in Fig. El and listed in section |31 We separately calculated the degree 
of coherence of each GRN over the full set of fixed points/cycles (a c ) and over 
the biologically relevant subset ( a \”°), both adopted from respective references. We 
compared them with the mean (< cC c and ) obtained from the associated ensemble of random 
networks described in Section 12.21 


Organism/Cell type 

a c (uniformly weighted) 

ol c (basin-size weighted) 


OL c 

<° 

^rand 

Ct c 

OL c 

or 

■ rand 

lx c 

S. pombe 

0.95 

1 

0.94 ± 0.04 

0.99 

1 

0.93± 0.05 

Mammalian cell 

0.88 

0.88 

0.85 ± 0.07 

0.88 

0.88 

0.86 ± 0.07 

S. cerevisiae 

0.97 

1 

0.80 ± 0.08 

0.99 

1 

0.75± 0.10 

A. thaliana whorls 

1 

1 

0.97 ± 0.04 

1 

1 

0.96 ± 0.04 

Myeloid progenitor 

0.89 

0.89 

0.82 ± 0.09 

0.90 

0.88 

0.84 ± 0.09 

Th-lymphocyte 

0.99 

0.96 

0.94 ± 0.05 

0.94 

0.92 

0.94 ± 0.05 


Table 2: Degrees of coherence associated with the regulatory networks shown in Fig. El 
af , o%°, a r c and , denote the coherence value calculated by using Eq.(JHl) over all attractors 
of the model dynamics, only the attractors observed in the wild type, and all networks 
in the corresponding randomized ensemble, respectively. 


Our central observation is that, despite the variability in k and p, all biological 
networks are more coherent than random networks constructed with similar structural 
parameters. Fig. El also shows the distribution of a c over the random ensembles for a 
better judgment of the bias towards coherence. The difference in coherence between the 
biological network and the random ensemble is visible but within acceptable bounds for 
most isolated examples. However, the fact that all display a positive deviation from 
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Figure 3: Degree of coherent regulation for biological vs random networks. The 
histograms of a c for random ensembles of 10 4 networks having similar structure is shown 
with their mean values (green line). The red line indicates the a c of the biological 
networks, whose exact values can be seen on Table El 


the respective ensemble means suggests an overall preference towards coherence. We 
quantified the statistically significance of this bias by using Fisher’s method [59 J on 
the present data, which yields a likelihood of 0.044 and 0.075 for a chance encounter 
of a this much or larger deviation in uniformly and basin-weighted cases, respectively. 
The degree of selection is appreciable, considering that the compared ensemble networks 
have not only the same number of nodes, the same number of interactions, and the same 
proportion of repressory/activatory regulation globally, but also the same local pattern 
of incoming and outgoing edges for each node. 

5. Coherence vs structural/dynamical network properties 

The bias observed in biological regulatory networks above provides motivation to 
investigate the relationship between coherence and other architectural or dynamical 
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Figure 4: Dependence of coherence value (a c ) on the network parameters (A:, p) with 
n = 11. As k increases from 1 to 10, the p value where a c is minumum for the fixed k 
shifts from 0.72 to 0.6. Coherence values corresponding to the biological networks are 
shown by spheres. 


determinants of a network. Identifying such connections could help one recognize 
coherent systems from certain telltale patterns instead of requiring detailed information 
about their function, and/or design them by means of simple guiding principles. 

Edge number & type: 

We first investigated the dependence of average coherence a c on the number of 
incoming/out going edges per node, A:, and the fraction of activatory interactions, p. 
To this end, we simulated the majority-rule dynamics in Ref. m on different (k,p) 
pairs equally spaced in the rectangle 0.909 < A: < 10 0 0 < p < 1, with n = 11. 
The result is shown in Fig. HI Each data point in Fig. H]is an average over 10 3 random 
networks which are generated by random shuffling of edges. No constraint was imposed 
on the shuffling process, apart from the requirement of connectedness (every node is 
accessible from every other node in the undirected network). 

One observes that, independent of the average connectivity k , random networks 
display minimum coherence when approximately two-thirds (60 — 72%) of the 
interactions are activatory. An intuitive understanding of this behavior was proposed 
in Ref. [30J where a subset of the results in Fig. H] (for a single k value) was reported. 

Coordinates corresponding to the studied biological networks are shown by colored 
spheres on the same figure. It is interesting that all of them are situated on the low-p 
slope of the minimum coherence valley. 

Number of active genes: 

We next asked if the number of active genes ( n a ) at a fixed point/cycle is a determining 
factor for coherence. After all, coherence reflects the harmony between the regulatory 















Coherent organization in gene regulation: a study on six networks 


11 


messages originating from these genes. To this end, we considered the attractors found 
within the ensembles generated by shuffling the edges of each GRN in Fig. El We 
grouped them according to n a and calculated the mean coherence within each group. 
The results in Fig. El show that attractors with higher number of active nodes are 
generally less coherent. This behavior may be understood by the intuitive fact that it 
is more difficult to reach consensus in a large group than in a small one. We remind 
that, this relationship is valid only within the restricted ensembles generated by the 
shuffling process described in Section El One can not speak of a monotonic dependence 
of coherence on n a if, for example, the difference in n a between two networks is due to 
different p values. This is evident from the fact that h a (p) calculated over the network 
ensembles used in Fig. 0] trivially increases with p while a c (p) is nonmonotonic for any 
fixed k. 



0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Coherence (a c ) 


Figure 5: Average ratio of active nodes in the point-attractors vs. average coherence 
in ensembles of model networks. The average is calculated over 10 4 networks chosen 
randomly from the null-model ensemble described in section 12.21 


Basin size: 

It is implicit from Table El and Fig. [3] that, coherent attractors are not always 
those with larger basins. We explicitly examined the variation in coherence as a 
function of the relative basin size of the attractors. The results (Fig. |6j), obtained 
over the networks in the randomized ensembles of each GRN separately, suggest no 
consistent relation between the basin size and coherence. Likewise, biologically relevant 
attractors [38 ; 1391 1371 , 55, [56, EZJ of GRNs in Fig. El are not always those with the largest 
basins. This is hardly surprising, since the initial state prior to differentiation or cell 
division is never determined randomly. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 


Basin size 

Figure 6: Average coherence vs. average basin size for the attractors (point or cycles) 
of 10 4 networks with similar structure. Note that, proposed rules of dynamics in some 
models appear to prohibit appearance of excessively dominant attractors. 


6. Conclusion 

We investigated the degree of coherent regulation in several regulatory networks 
responsible for cell cycle and differentiation in various organisms or cell types, by 
means of a recently proposed measure of coherence. We found that, even though most 
networks are moderately more coherent than expected (in reference to architecturally 
similar network ensembles), cumulatively there is a statistically meaningful bias towards 
coherence. Our findings lend support to the thesis that pressure for coherent regulation 
is one of the factors driving the evolution of GRNs. It is not difficult to imagine 
a Hebbian-like mechanism |60j in the given context, where regulatory interactions 
incompatible with the required expression profile get eliminated over time. The fact that 
transcriptional regulatory networks have been observed to have the fastest mutation rate 
among various biological networks |61j suggests that their evolutionary dynamics may 
be sufficiently fluid for effective Hebbian selection. 

We also showed, on randomly generated model GRNs, that coherence is harder 
to achieve with increasing number of interactions per node, and with an inhibitory 
interaction ratio around 1/3. Furthermore, coherent networks typically involve a smaller 
number of active genes at their steady states, compared to arbitrary networks with 
similar edge composition and local connectivity. It could be interesting to focus on 
networks following the identical expression pattern in time as, say, the yeast cell-cycle 
(see, e.g., Ref. |62j) and check if those with higher coherence are more similar to their 
biological counterpart or encompass other desirable properties such as higher robustness 
or better controllability. Finally, from an architectural perspective, coherence can be 
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viewed as a design choice for any natural or artificial network where inhibitory and 
activatory interactions coexist. One might then ask how to build a coherent system from 
scratch, or how to enhance coherence in an existing system with minimal intervention. 
We hope our findings to trigger further theoretical investigations in these directions. 
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Below, we list the rules of regulatory dynamics for each GRN model adopted from the 
literature. 


Network 

Gene 

Boolean Update Function 



GATA-2 

GATA-2 A (GATA-1 A FOG-1) A PU.l 



GATA-1 

(GATA-1 V GATA-2 V Fli-1) A PU.l 



FOG-1 

GATA-1 



EKLF 

GATA-1 A Fli-1 


Myeloid 

Fli-1 

GATA-1 A EKLF 


differentia- 

SCL 

GATA-1 A PU.l 


tion 

C/EBP a 

C/EBP q A (GATA-1 A FOG-1 A SCL) 



PU.l 

(C/EBP a V PU.l) A (GATA-1 V GATA-2) 



cJun 

PU.l A Gfi-1 



EgrNab 

(PU.l A cJun) A Gfi-1 



Gh-1 

C/EBP a A EgrNab 



Table 3 

Update functions for myeloid differentiation. 

Network 

Gene 

Boolean Update Functiont 


CycD 

CycD 


Rb 

(CycD A CycE A CycA A CycB) V (p27 A CycD A CycB ) 


E2F 

(Rb A CycA A CycB) V (p27 A Rb A CycB) 


CycE 

E2F ARb 

Mammalian 

CycA 

(E2F ARb A Cdc20 A (Cdhl A Ubc)) 


cell cycle 

V (CycA A Rb A Cdc20 A (Cdhl A Ubc)) 



p27 

(CycD A CycE A CycA A CycB) V ( p27 A (CycE A CycA ) A CycB A CycD ) 


Cdc20 a 

CycB 


Cdhl 

(CycA A CycB) V Cdc20 V (p27 A CycB) 


UbcHIO 

Cdhl V (Cdhl A Ubc A (Cdc20 V CycA V CycB)) 


CycB 

Cdc20 A Cdhl 


Table 4: Update functions for mammalian cell cycle. 
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Update Function^ 


Network 


Node i 


Gene 


S. cerevisiae 
cell cycle 


1 

2 

3 

4 

5 


Cln3 

Clnl,2 

Cdc20 

Mem 

Swi5 


Si(t + 1) — 



6 

7 

8 

9 

10 
11 


SBF 

MBF 

Sicl 

Clb5 

Cdhl 

Clbl 


Si(t + 1) — 


1, / > 0 
0, I < 0 

Si(t), 1 = 0 


S. pombe 
cell cycle 


1 

2 

3 


SK 

SLP 

PP 


Si(t + 1 ) — 



4 

5 

6 

7 


Ste9 

Ruml 

Cdc2 

Cdc2* 

Weel 

Cdc25 


Si(t + 1) — 


1 , I > 0 t 

0 , I <9i 

Si(t ), l = 0i, 


A. thaliana 
whorl 

differentiation 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


EMF1 

TFL1 

LFY 

API 

CAL 

LUG 

UFO 

BFU 

AG 

AP3 

PI 

SUP 


Si(t +1) — 



In the update functions I = JU = i A l] S l {t) where A is the weighted adjacency matrix 
and 6 's indicate the thereshold values for the nodes. For detailed explanations for these 
values one can refer to EtEHlESI- _ 


Table 5: Update rules of the regulatory dynamics for the budding yeast ( Saccharomyces 
cerevisiae ), the fission yeast ( Saccharomyces pombe), and the flower Arabidopsis thaliana. 
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Network 

Gene 

Boolean Update Function^ 

Lymphocyte 

Differentiation 

C^ 

i 

Z 

Pm 

i— i 

^(9), K^ll), AT X (9,11) 

IL-4 

^2(12) 

IL-12 


IFN-yR 

A' 4 (l) 

IL-4R 

K 5 (2) 

IL-12R 

M3) 

STAT1 

K,( 4) 

STAT6 

K 8 (5) 

STAT4 

K g (6) 

SOCS1 

K 10 {7),K 10 { 11), ^(7,11) 

T-bet 

Ku{7),Ku(U),K n (7,n) 

GATA-3 

K 12 (8) 

K, ( { j }) is the set of nodes j that when simultaneously active turn the node i on. For 
details one can refer to |H7|. 


Table 6: Update functions for lymphocyte differentiation. 
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